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It may be possible to derive a constituent approximation for bound states in 
QCD using hamiltonian light-front field theory. Cutoffs that violate explicit gauge 
invariance and Lorentz covariance must be employed. A similarity renormalization 
group and coupling coherence are used to compute the effective hamiltonian as an 
expansion in powers of the canonical QCD running coupling constant. At second 
order the QCD hamiltonian contains a confining interaction, which is being studied 
using bound state perturbation theory. Explicit constituent masses appear because 
of symmetry violations, and confinement also produces mass gaps, leading to the 
possibility of an accurate non-perturbative constituent approximation emerging in 



CL|. light- front QCD. 

< 

1. Introduction 

The solution of Quantum Chromodynamics in the non-perturbative domain remains 
' one of the most important and interesting unsolved problems in physics. QCD is be- 

lieved to be the fundamental theory of the strong interaction, but even its definition in 
the non-perturbative domain is problematic. There are many sources of difficulty, but 
they can all be traced to the fact that QCD is formulated as a theory of an infinite 
' number of degrees of freedom that span an infinite number of energy scales. 

\Q . The basic assumption upon which all of our work is based is that it is possible 

0^ ' to derive a constituent picture for hadrons from QCD [1-5]. If this is possible, non- 

r-j ■ perturbative bound state problems in QCD are approximated as coupled, few-body 

fS |' Schrodinger equations. For a meson, we then have, 

Q_|' p2 , M 2 

jg ■ P~ I *) = ± p + I *>. (1) 

where, 

*} = <f>qq I 99> + <f>qqg I 99.9} + ' ' (2) 

H , 

P is the light-front hamiltonian, P± is the total transverse momentum, P + is the 
total longitudinal momentum, and M is the invariant mass of the state. We assume 
that to 'leading order' a low-lying meson can be approximated as a quark/antiquark 
pair, with additional quarks and gluons producing 'perturbative' corrections that can 
be systematically computed. 

Many severe problems must be overcome to arrive at this formulation of the bound 
state problem; however, the final advantages are huge. The result is a formulation of the 
non-perturbative problem in a form directly accessible to physical intuition, which has 
proven essential for guiding approximations in atomic calculations. Variational methods 



1 



and large matrix diagonalization are powerful numerical tools that can be used after 
the hamiltonian is determined. 

I must emphasize that it is not our intent to simply force the constituent approxima- 
tion on the theory by employing a Tamm-Dancoff truncation on the number of particles 
ab initio. We worked on such an approach initially [6] , and gained valuable insights; but 
it became clear that we have no good method of controlling the nonlocalities resulting 
from particle number truncation without a dynamical mechanism that naturally limits 
the number of particles in a state. 

Any student of field theory should immediately be suspicious of the possibility that 
a constituent approximation can arise, although QED provides an important accepted 
example that guides much of our work [2] . How can a constituent approximation arise 
in any field theory? 

Fock space is extremely large, an infinite sum of cross products of infinite dimensional 
Hilbert spaces. It is not obvious that the low-lying eigenstates should have significant 
support only in the few-body sectors of Fock space. In fact, this simply does not 
happen in perturbation theory. In perturbation theory high-energy many-body states 
do not decouple from low-energy few-body states. Consider an electron mixing with 
high-energy electron/photon states. The error made by simply throwing away the high 
energy components of the state is infinite. Moreover, there are an infinite number of 
scales and both the electron and photon that 'dress' the low-energy bare electron are 
in turn dressed by additional pairs, ad infinitum. 

The lesson here is quite old. Without regularization and renormalization a con- 
stituent picture is impossible. Renormalization may allow us to move the dynamical 
effects of high-energy, many-body states from the eigenstate to effective interactions 
between effective quarks and gluons. 

Low-energy many-body states also do not decouple from low-energy few-body states. 
In fact, it is common lore that hadrons are excitations on an extremely complicated 
vacuum. Students of QCD expect the infinite-body vacuum to be an integral part of 
every hadron eigenstate. This is the problem that leads us to use light-front coordinates, 
just as it motivated the use of the infinite momentum frame for the formulation of 
the parton model. In light-front coordinates physical particle trajectories satisfy the 
kinematic relativistic constraint 

P + > , (3) 

because all velocities are equal to or less than the velocity of light. Since longitudi- 
nal momentum is conserved, the only states that can mix with the zero momentum 
bare vacuum are those in which every bare parton has identically zero longitudinal 
momentum. For a free particle of mass to, the light-front energy is 

„- = Ei±^. ( 4) 

This energy diverges as p + approaches zero, which must happen as the number of 
particles grows for fixed total longitudinal momentum. Thus, in light-front coordinates 
all many-body states become high energy states, leading us back to the original problem 
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of replacing the effects of high energy states with effective interactions. This argument 
is naive, but there is little profit in elaborating further at this point. 

Finally, manifest gauge invariance and manifest covariance apparently require all 
states to contain an infinite number of particles. This is most easily seen, for example, 
by considering rotation operators. Rotations are dynamical in light-front coordinates 
and the generators contain interactions that change particle number. No state with 
a finite number of particles transforms correctly under rotations. We use cutoffs that 
violate these symmetries, which must then be repaired by effective interactions that 
remove all cutoff dependence. The constituent approximation is possible only if these 
symmetries are also treated approximately. Proposing the violation of manifest gauge 
invariance is heresy in the QCD community, but heresy sometimes leads to progress in 
science. 

There is a long list of questions concerning how a constituent approximation can arise 
in QCD, but I mention only one. How can confinement emerge without a complicated 
vacuum? Since we have a hamiltonian, we can use a variational calculation to study 
what happens as a quark/antiquark pair are separated to infinity. Since the addition 
of gluons can only lower the energy, we must find that 

(<Pqq I H | <j> q q) > OO aS R ► OO . (5) 

Here R is the quark separation, and the only way this matrix element can diverge is if 
the hamiltonian contains a two-body interaction that diverges. We will see below that 
this apparently happens. 

This discussion is not intended to convince the skeptical reader that a constituent 
approximation is valid. However, the assumption that a constituent picture emerges 
from QCD provides strong guidance. A hamiltonian approach is indicated. Cutoffs that 
limit mixing of high and low energy states are required, and they must violate explicit 
rotational covariance and gauge invariance. All non-perturbativc effects attributed to 
the vacuum in other approaches must directly appear in few-body effective interactions. 

2. Light-Front Renormalization Group 

The renormalization of the hamiltonian and all other dynamical observables begins 
with the observation that no physical result can depend on the cutoff. In the Schrodingcr 
equation, 

p 2 _i_ M 2 

Pa I *a> = ± p + ' I *a) , (6) 

the eigenvalue, M, cannot depend on the cutoff. The hamiltonian, P~, must depend 
on the cutoff, as does the eigenstate. Wilson's renormalization group was formulated 
starting with the observation that physical matrix elements cannot depend on the cutoff, 
and we have adapted his approach to the light-front problems we face [7] . 

It is not possible to discuss cutoff-independence if the cutoff is fixed, so the central 
operator in Wilson's renormalization group is a transformation that lowers the cutoff. 
Given a transformation, T, that lowers the cutoff by a factor of 1/2, for example, we can 
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define a renormalized hamiltonian to be one which has a finite cutoff but results from 
an infinite number of transformations. The transformation determines what operators 
must be precisely controlled for this limit to exist. Near a fixed point (i.e., a hamiltonian 
that does not change under the transformation), these operators can be classified as 
relevant and marginal. 

In the perturbative regime relevant and marginal operators are determined by their 
naive engineering dimension. In light-front field theory there is no longitudinal locality, 
only transverse locality, so it is the transverse dimension of an operator that determines 
its classification. However, while there are a finite number of relevant and marginal 
operators in equal-time field theory, the violation of longitudinal locality in light-front 
field theory implies that ratios of longitudinal momenta can appear, allowing entire 
functions of longitudinal momentum fractions to appear in each relevant and marginal 
operator. At first sight this appears to be a disaster; however, one paradox we faced 
above was how complicated interactions associated with non-perturbative effects such as 
confinement could arise in few-body operators. This is possible because of the violation 
of longitudinal locality. 

To develop a light-front renormalization group we must decide what cutoff to im- 
plement and then derive a transformation that runs this cutoff. It is possible to use a 
cutoff on the total invariant-mass of states, as is commonly done in DLCQ calculations 
for example; however, such cutoffs lead to strong spectator dependence and small en- 
ergy denominators appear in the resultant effective interactions. We use a cutoff on the 
change in free energy. If the hamiltonian is viewed as a matrix such a cutoff limits how 
far off diagonal matrix elements can appear. 

There is not enough space to elaborate the transformation that runs this cutoff, so 
I refer the reader to the literature [1,2,4,8,9]. The transformation is unitary, leading to 
what Glazek and Wilson call a similarity renormalization group [8,9]. Let H = ho + v, 
where h is a free hamiltonian with ho \ fa) = E oi \ fa), and v is cut off so that 

(fa\v\fa) = , (7) 

if \E 0i — E j \ > A. If this cutoff is lowered to A', the new hamiltonian matrix elements 
to 0{v 2 ) are 



H 'ab = (4>a\h0 + v\(f> b ) -^VakVkb 



0(|A nfc |-A')0(|A nfc |-|AM|) 

Eok — Eoa 

0(|A 6fe |-A')0(|A 6fe |-|A afe |) 



+ 



Eok — Eob 



(8) 



where A^ = Eoi — Eqj and \Eo a — Eob\ < A'. To follow the details of the discussion it 
is important to remember that there are implicit cutoffs in this expression because the 
matrix elements of v have already been cut off so that Vij = if \E ai — E j \ > A. 

It is rather easy to understand this result qualitatively. We have removed the cou- 
pling between degrees of freedom whose free energy difference is between A' and A, 
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so the effects of these couplings are forced to appear in the new hamiltonian as direct 
interactions. To first order, the new hamiltonian is the same as the old hamiltonian, 
except that couplings of states with energy differences between A' and A are now zero. 
To second order, the new hamiltonian contains a new interaction which sums over the 
second-order effects of couplings that have been removed. The second-order term in the 
new hamiltonian resembles the expression found in second-order perturbation theory, 
which is not surprising since the new hamiltonian must produce the same perturbative 
expansion for eigenvalues, cross sections, etc. as the original hamiltonian. 

Equation (9) shows how the hamiltonian changes when the cutoff is lowered, and 
the next step is to determine from this change what hamiltonians can emerge from an 
infinite number of transformations. The simplest result is a fixed point hamiltonian, one 
which does not change under the transformation. In 3 + 1 dimensions the only known 
fixed points are free field theories. Coupling coherence is a generalization of the fixed 
point idea [2,7,10]. A coupling coherent hamiltonian reproduces itself in form, but one 
or more couplings run while all additional couplings are invariant functions of these 
running couplings. For example, in QCD the canonical coupling runs at third order. To 
second order in this coupling, all interactions must reproduce themselves exactly, with 
A — > A'. It is not trivial to implement this simple-sounding constraint, but at each order 
it determines the hamiltonian. In all calculations to date the resultant hamiltonian is 
unique, and all broken symmetries are restored to the order at which the hamiltonian 
is fixed. 

To second order a generic coupling coherent hamiltonian that contains v must also 
contain, 



Hab = (<Pa\h + v\<p b ) -^2v ak Vkb 



fl(|A afc |-A)0(|A afc |-|A 6fc |) 
9(\A bk \-A)9(\A bk \-\A ak \) 



+ 



Eok — Eob 



(9) 



or 



H a b — 



b a \h + v\(j) b ) +^2v ak v kb 



e(A-|A ofc |)fl(|A ofc |-|A Wb |) 



+ 



0(A 



Eo k — Eoa 

|A 6fc |)0(|A 6fc | 



|A afe |) 



Eok - E ob 



(10) 



Note that v in these expressions is the same as that above only to first order. The 
coupling coherent interaction in H is written as a power series in v which reproduces 
itself under the transformation, except the cutoff changes. In higher orders the canonical 
variables also run. 

The light-front similarity renormalization group and coupling coherence fix the QCD 
hamiltonian as an expansion in powers of the running canonical coupling. 
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3. QCD: A Strategy for Bound State Calculations and Confinement 

While realistic calculations will no doubt require a more elaborate procedure, a 
relatively simple strategy for doing bound state calculations can now be outlined [2,4]. 

i) Start with the canonical hamiltonian, H can , and use the similarity renor- 
malization group and coupling coherence to compute, 

H A = h A + g A h A + g 2 A h A + --- (11) 
Truncate this series at a fixed order. 

ii) Choose an approximate hamiltonian that can be treated non- 
perturbatively, 

H A = H A + V A . (12) 
You must choose A and Hq to minimize errors. 

iii) Accurately solve Hq as the leading approximation. 

iv) Compute higher order corrections from V A using bound state pertur- 
bation theory. 

v) To improve the calculation further return to step (i) and compute the 
hamiltonian to higher order. 

There are two principal reasons that this strategy will fail for QCD. First, the 
hamiltonian is computed perturbatively so that errors in the strengths of all operators 
are at least as large as a power of a. Small errors in the strengths of irrelevant operators 
tend to produce even smaller errors in results. However, errors in marginal operators 
tend to produce errors of the same order in results and small errors in relevant operators 
tend to produce exponentially large errors in results. At the minimum we expect that 
we will have to fine tune relevant operators, which means tuning a finite number of 
functions of longitudinal momenta. Second, chiral symmetry breaking operators (where 
light- front chiral symmetry should be distinguished from equal-time chiral symmetry [1]) 
will not arise at any order in an expansion in powers of the strong coupling constant. Wc 
must work in the broken symmetry phase of QCD ab initio and insert chiral symmetry 
breaking operators. Simple arguments lead us to expect that only relevant operators 
need to be considered if transverse locality is maintained, but there are no strong 
arguments for transverse locality in these operators. 

Despite these limitations, this strategy may be applied to the study of bound states 
containing at least one heavy quark [5], as discussed by Martina Brisudova in these 
proceedings; although even here masses should be tuned, as expected. The strategy is 
conceptually simple and there are no ad hoc assumptions. 

The first step is to compute the effective QCD hamiltonian to order a. I refer the 
reader to the literature for details on the canonical hamiltonian [11]. The first applica- 
tions of the approach are to mesons [5] , and we assume that for sufficiently small cutoffs 
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we can choose Hq to contain only interactions in H A that do not involve particle pro- 
duction or annihilation, as dictated by our initial assumption that a constituent picture 
will arise. This means we can first focus on operators that act in the quark/antiquark 
sector. I emphasize that all operators must be computed, and without the confining 
interactions in sectors containing gluons the entire approach would make no sense. 

First consider the second-order correction to the quark self-energy. This results 
from the quark mixing with quark-gluon states whose energy is above the cutoff. If we 
assume that the light-front energy transfer through the quark-gluon vertex must be less 
than A 2 /V + , the coupling coherent self-energy for quarks with zero current mass is 



Here the quark has longitudinal momentum p + , while the longitudinal momentum scale 
in the cutoff is V + . The first and most interesting feature of this result is that I have 
been forced to introduce a second cutoff, 



which restricts how small the longitudinal momenta of any particle can become. With- 
out this second cutoff on the loop momenta, the self-energy is infinite, even with the 
vertex cutoff. This second cutoff should be thought of as a longitudinal resolution. As 
we let e — > we resolve more and more wee partons, and in the process we should 
confront effects normally ascribed to the vacuum. In this case the wee gluons are re- 
sponsible for giving the quark a mass that is literally infinite. Theorists who insist on 
deriving intuition from manifestly gauge invariant calculations may find this interpre- 
tation repugnant, but within the framework of a light-front hamiltonian calculation it 
is quite natural. 

This second, infrared cutoff poses a problem. If we introduce a second cutoff, 
shouldn't we introduce a second renormalization group transformation to run this cut- 
off and find the new counterterms required by it? No. / will insist that all divergences 
associated with e — > cancel exactly in all physical results for color singlet states. The 
important question is how can these divergences cancel so that mesons have a finite 
mass, and the answer to this question leads to confinement. 

A nearly identical calculation leads to the second-order self-energy of gluons, and 
the dominant term goes like 



The quark and gluon masses are infinite, which is half of the confinement mechanism. 
In addition to one-body operators we find quark-quark, quark-gluon, and gluon-gluon 
interactions. As we lower the cutoff, we remove gluon exchange interactions, and these 
are replaced by direct interactions. The analysis of all of these interactions is nearly 




(13) 



(14) 




(15) 
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identical, and I consider only the quark- antiquark interaction. This interaction includes 
two pieces, instantaneous gluon exchange which is in the canonical hamiltonian, and an 
effective interaction resulting from high energy gluon exchange. To study confinement 
we need to examine the longest range part of the total interaction, which is a piece that 
diverges as longitudinal momentum exchange goes to zero. I outline the calculation 
[2,4]. 

High energy gluon exchange cancels part of the instantaneous gluon exchange inter- 
action, leaving 



V singular = -4g 2 A C F ^/p+p+k+k+ 9(A 2 /V + -\q-\) 9(\q+\-eV + ).(16) 

Here the initial and final quark (antiquark) momenta are p\ and P2 (fci and £2), and 
the exchanged gluon momentum is q. The energies are all determined by the momenta, 
p^ = p\ x jp\, etc. This part of the interaction is independent of the spins. If A « 
A-qcd, we expect further gluon exchange to be suppressed, and we are left with this 
singular interaction between the quark and antiquark. 

The next step in the analysis is to take the expectation value of this interaction 
between arbitrary quark-antiquark states, (vf^ I V | ^1). If we define 

Q = — g — ' 9 = Pl~P2, (17) 

and expand the wave functions about q = 0, we find a divergence in the expectation 
value, 



= -"i^r^O J dQ+ JJ^ 2(Q)MQ)- (is) 

Unless 0i and 4>2 are the same, this vanishes by orthogonality. If they are the same, 
this is exactly the same expression we obtain for the expectation value of the quark plus 
antiquark divergent mass operators; except with the opposite sign. Therefore, there is 
a divergence in the quark-antiquark interaction that is independent of their relative 
motion and which exactly cancels the divergent masses! These cancellations only occur 
for color singlets, and they occur for any color singlet state with an arbitrary number 
of quarks and gluons. Moreover, these cancellations appear directly in the hamiltonian 
matrix elements, so we can take the e — > limit before diagonalizing the matrix. 

This is half of the simple confinement mechanism. At this point it is possible to 
obtain finite mass hadrons even though the parton masses diverge. However, since 
the cancellations are independent of the relative parton motion, we must study the 
residual interactions to see if they are confining. Since I am interested in the long-range 
interaction, I will study the fouricr transform of the potential and compute V(r) — V(0) 
so that the divergent constant in which we are no longer interested is canceled. 

2 n ^2 

V singu i ar (r) ~ V singu i ar (0) -> A ^ log(|a;~|) , (19) 
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when x± — and \x \ 



oo; and 




g\c F k 2 



iog(M) , 



(20) 



2n 2 P+ 



when — > cxo and x~ = 0. This potential is not rotationally symmetric, but it 
diverges logarithmically in all directions. 

If the potential is not rotationally symmetric, how can rotational symmetry be re- 
stored? In light-front field theory rotations are dynamical. While it may be possible for 
rotational symmetry to be realized approximately in low-lying quark-antiquark states, 
exact symmetry requires additional explicit partons and even approximate rotational 
symmetry will require additional partons if we study highly excited states. We expect 
excited physical states in which a quark and antiquark are separated by a large dis- 
tance to contain gluons. There is no reason to assume that the gluon content of these 
states is the same when the state is rotated, so rotational symmetry will be restored in 
highly excited states only if we allow additional partons. This complicates our attempt 
to derive a constituent picture, but we only need the constituent picture to work well 
for low-lying states. The intermediate range part of the potential is rotationally sym- 
metric, and we may expect the ground state hadrons to be dominated by the valence 
configuration. 

Isn't the confining potential supposed to be linear and not logarithmic? There is 
no conclusive evidence that the long-range potential is linear, and heavy quark phe- 
nomenology shows that a logarithmic potential can work quite well; however, lattice 
calculations provide strong evidence for a linear potential. However, low-lying states 
are not sensitive to the longest-range part of the interaction; and light quark-antiquark 
pairs prevent even exited states from being sensitive to the longest-range part of the 
interaction. In any case, I do not want to argue that these calculations show that the 
long-range potential in light-front QCD is logarithmic. Higher order corrections could 
produce powers of logarithms that add up to produce a linear potential. 

The above argument seems to apply directly to QED at first sight. Is QED confining? 
There is a confining interaction between charged particles in the hamiltonian, but there 
is no strong interaction between charged particles and photons. To see if confinement 
survives in QED we should include the confining interaction in Ho, and then compute 
corrections in bound state perturbation theory. In QED the second order correction 
from photon exchange below the cutoff exactly cancels confinement. This implies that 
if confinement is included in Ho, higher order corrections are large. If the Coulomb 
interaction, which appears in H also, is included in H , bound state perturbation theory 
appears to converge rapidly. On the other hand, in QCD gluons also experience a 
confining interaction. When second order bound state perturbation theory is used to 
study the effect of the exchange of confined gluons, it is seen that gluon exchange does 
not cancel the confining interaction in Ho; so this picture of confinement is at least 
self-consistent. 

The important point is that H contains a confining interaction that we are free to 
include in Ho, giving us some hope of finding a reasonable bound state perturbation 
theory for hadrons that resembles the bound state perturbation theory that has been 
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successfully applied to the study of atoms. 
Summary 

A constituent picture of hadrons may emerge in QCD if we use: 

• hamiltonian light-front field theory 

• a cutoff of order Aqcd on energy changes that violates manifest covariance and 
gauge invariance 

• a similarity renormalization group and coupling coherence 

Bound states can be studied using bound state perturbation theory. The effective 
hamiltonian is computed as an expansion in the strong coupling, and then divided into 
H = Hq + V, with V treated perturbatively. Hq must include all essential interactions. 
We have found that H contains an order a logarithmically confining two-body interac- 
tion between all colored partons, and we have begun studies of bound states using this 
confining interaction, as discussed in the talk by Martina Brisudova. 
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